The microscopic meaning of grand potential resulting from combinatorial approach to 

a general system of particles 
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We present a completely new approach to the problem of interacting fluids, which we believe may 
provide important insights into microscopic mechanisms that lead to the occurrence of phase transi- 
tions. The approach exploits enumerative properties and combinatorial meaning of Bell polynomials. 
We derive the exact formula for probability of a general system of N particles at temperature T 
to consist of k weakly coupled clusters of various sizes. We also show that the grand potential of 
the system may be considered as the exponential generating function for the number of internal 
states (thermodynamic probability) of these clusters. The microscopic interpretation of the grand 
potential is novel and surprising, especially if one recalls that until now the only thermodynamic 
meaning of this free energy was known. We also derive an approximated expression for the density 
of states. 

PACS numbers: 05.20.-y, 05.20.Gg, 02.10.Ox 



I. INTRODUCTION AND MOTIVATION 

This work springs from a fascination with the con- 
cepts of equilibrium, interaction, and collective behavior 
in physical systems. The analyzed systems are consid- 
ered to be made up of separate clusters which are com- 
posed of particles (elements, units etc.) interacting via 
short-range potentials. Exemplary systems where the as- 
sumptions apply are imperfect gases, for which the total 
interparticle potential energy is given by the sum of pair 
interactions characterized by a hard core and a weak at- 
tractive region. In this study we assume a particular 
type of interaction between clusters termed weak cou- 
pling. Weak coupling means that the only interaction 
between clusters is that provided by the migration of par- 
ticles between them. For this reason, within the frame- 
work of statistical mechanics, the system can be actu- 
ally described as consisting of disjoint, non-interacting, 
and independent clusters. A systematic method of deal- 
ing with such systems was first developed by Ursell 
whose work was later extended by Mayer [2t3|; due to 
whom are the major advances in this field. 

With this paper we start a new line of theoretical re- 
search on imperfect gases (or interacting fluids) . Our ap- 
proach originates in basic concepts of statistical physics 
(ensemble theory) [8| and enumerative combinatorics [9j . 
In particular, we derive a new, exact and general combi- 
natorial formula for probability of a system of N particles 
at temperature T to consist of k clusters of various sizes. 
To achieve the result we use the generating function for 
Bell polynomials [Tfjj . To interpret it we exploit combina- 
torial meaning of this function as given by the so-called 
exponential formula fllT ] , which is the cornerstone of the 
art of counting in enumerative combinatorics. 

To be definite, the exponential formula deals with the 
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question of counting composite structures that are built 
out of a given set of building blocks [ll[ (in the exemplary 
case of imperfect gases, the term "composite structure" 
stands for "microstate" of the system, while the term 
"block" means "cluster"). If you know how many blocks 
of each size there are, the formula tells you how many 
composite structures of each size can be built from those 
blocks. In short, the exponential formula states that the 
exponential generating function, F(x) = J2^=q fnX n /n\, 
for composite structures is the exponential of that for 
building blocks, f(x), i.e. F{x) = exp[/(x)]. With this 
statement in mind, the grand partition function in statis- 
tical physics [26j], H(/8, z) — exp[— /3<J>(/3, z)], which is also 
a kind of generating function for microscopic arrange- 
ments of the system investigated, gains a completely 
new interpretation as the generating function for com- 
posite states having a well-defined internal (clustered) 
structure. Furthermore, properties of this internal struc- 
ture are encoded in another generating function, $(/?, z), 
which appears to be free energy of the system. 

The outlined understanding of the grand partition 
function as given by the exponential formula is novel, but 
its theoretical basis is partially covered by the cluster ex- 
pansion method, set out over 50 years ago in the already 
mentioned works of Mayer I5|, Jfj . Mayer's ideas were 
later developed by others [ilj.Tl3l|7for example by the dis- 
tinguished physicist Uhlenbeck and his student Riddel, 
in whose thesis the exponential formula first appeared 
14]. However, following initial interest among physi- 
cists, the idea of enumerative uses of generating functions 
was generalized and extended mainly by mathematicians 
(cf . [1, [ll[ ) , and the physical origins of the concept have 
almost been forgotten. 

This happens because the Mayer's original approach 
to imperfect gases failed to match real systems. Fore- 
most it failed in the description of condensation effects 
in dense interacting gases for which it was originally de- 
veloped. Thus, is it worth to coming back to this un- 
successful story? We believe, it is wort it, especially 
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that given the present comprehension of different enu- 
merative methods, one can easy trace the Mayer's failure 
and hopefully explain the controversy over rather sub- 
tle mathematical points involved in his condensation and 
critical point theory. It would be the great success of 
equilibrium statistical mechanics, because although the 
alternative formalism of Yang and Lee fl5| makes pos- 
sible a simple mathematical description of condensation 
phenomena but it has still (after more than six decades 
of research) the disadvantage of being less detailed and 
physically rather obscure (see [7j, p. 123) in comparison 
with the Mayer's approach. 

Despite of its strong connections with the Mayer's the- 
ory, the ideas described in this paper were developed en- 
tirely independently. For the reason, the connections are 
not easily seen at first glance. In fact, for a long time of 
struggling with the exponential formula, Bell polynomi- 
als etc., we were unaware of the previous contributions 
in the field. Now, after several months of rummaging 
through literature on the topic, we can locate our studies 
as much more general than the Mayer's approach, and 
as being complementary to Markov models of systems in 
stochastic equilibrium studied by Whittle [16| , and nowa- 
days independently developed within the framework of 
stochastic combinatorial processes (l7l Il8j. 

The work reported in this paper is quite general and 
fairly abstract. Our main result is that the grand po- 
tential may be considered as the exponential generating 
function for the number of internal states of clusters of 
different sizes. It puts the grand potential (also called 
the Landau free energy), $(/3, z), in a completely new 
perspective. Up to now the only macroscopic interpre- 
tation of this thermodynamic state function was known 
and considered (see Sec. 2F5 in Q). The ideas advanced 
in this paper (and the forthcoming ones) reveal its mi- 
croscopic meaning. This, in turn, may shed new light 
on microscopic origins of phase transitions, since they 
are traditionally classified according to the lowest-order 
derivative of the thermodynamic potential which changes 
discontinuously (or diverges) at the transition point (see 
Chaps. 3 and 7 in |8j). In subsequent papers we will ap- 
ply the methods here outlined to lattice gasses. We will 
also show how the virial expansion 0, [l!|of the equation 
of state for imperfect gases can be simply derived from 
our approach. 

This paper is organized as follows. In the next sec- 
tion, Sec. [TT1 the main result of this paper, i.e. the ex- 
act combinatorial formula is derived for thermodynamic 
probability of a general system of N interacting parti- 
cles to consists of k weakly coupled clusters. Then, in 
the subsequent section, Sec. IIII1 properties of the for- 
mula are discussed, the general integral expression for 
the density of states is introduced, the hypothesis about 
microscopic meaning of the grand potential is advanced, 
and the example of the ideal gas is worked out. Some 
general conclusions are given in the end of the paper. 



II. DERIVATION OF THE MAIN RESULT 

In this paper, we limit ourselves to the simplest pos- 
sible case: a system composed of identical monoatomic 
molecules without internal degrees of freedom. The ther- 
modynamical state of the system is given by the tem- 
perature, T, and the chemical potential per molecule, fx. 
We assume that the classical treatment is adequate and 
use the formalism of the grand canonical ensemble (see 
Sec. 7G in [8() to describe the open system in the infinite 
volume limit. 

To derive the main result, we write the grand par- 
tition function (GPF) of a general system of particles 
in two ways. Firstly, we use a slightly modified text- 
book approach which exploits the definition of DOS as 
the number of states whose energies lie in some narrow 
range, between E and E + dE. Secondly, we make use of 
the generating function for partial Bell polynomials [To| . 
which has a very nice combinatorial interpretation given 
by the exponential formula [ll[ . Then, we compare both 
expressions. For the reasons of clarity GPF is denoted, 
respectively, by Si and S2, in the two approaches. 

In the infinite volume limit the grand partition func- 
tion, Si(/3, z), of a general system of particles is defined 
as the sum of Boltzmann exponential factors, e^( £ ~'' N ' , 
over all microscopic realizations (microstates), i7, of the 
system 

00 

Si 03, z) = J2 e-WM-""™ = 1 +J2 zNz (^ N), 

_ (1) 

where j5 = (fcsT) 1 , [i, E and N, have their usual mean- 
ings [27| . the parameter z = is called fugacity (or 
activity), and 

poo 

Z((3, N) = e -P E ( n ) — / g (E,N)e- fiE dE (2) 

stands for the canonical partition function (PF). Insert- 
ing Eq. ([2]) into Eq. (fTJ) and then multiplying the resulting 
expression by the sum over the probability distribution, 
f(k,E), that the considered system having energy E con- 
sists of k weakly coupled (i.e. statistically independent 
disjoint groups of particles) 

JV 

^f(k,E) = h(E) = l, (3) 
fc=i 

one gets 

3x09,*) = 

00 „oo / N \ 

= l+}2z N g(E,N)e-f ,E [J2f(k,E))dE (4) 

N=l J ° \k=l ) 

00 poo N_ 

= l+J2* N Y,f^ E ^9(E,N)e^ E dE. (5) 

N=l •'° k=l 
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To proceed further we have to assume that [20|: i. the 
sum in Eq. © is identically equal to one almost every- 
where, that is, except on a set of measure zero, and ii. 
the unknown function, f(k, E), multiplied by probability, 
g(E, N)e~ l3E , that energy of the considered system of N 
particles lies in a narrow range between E and E + dE, 
is an essentially bounded function of E. In line with the 
two assumptions, the grand partition function, Si, can 
be rewritten as 

oo N .qo 

N=l 

(6) 

The grand partition function, S2(/3, z), may also be 
written in the following way 



-/8*09,z) 



exp 



z 

^-^ ml 



(7) 



where <f>(/3,z) is the grand thermodynamic potential, 
derivatives of which, <p m {j3), are denoted by 



,00) 



with 



<h(P) = 0, 



(8) 



(9) 



which results from simple comparison of Eqs. (p} and 
for z — 0. Now, taking advantage of the generating 
function for partial Bell polynomials, -B/v.fc({</>n}), which 
are the polynomials in an infinite number of variables 
{</>„} = </>!, 09, </>3 ... , defined by the formal double series 
expansion [lOj 



G{t,u) 



exp 



(10) 



V 4>m— r 
t-^ m! 

m— 1 
N=l ' fc=l 

Eq. (JTJ) can be rewritten as 

oo JV 

3 2 (/3, z) = 1 + £ — "£(-P) k B N , k ({<t>nWM (12) 
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where B Njk ({<j) n }) = B N , k {4>i,4>2, ■ • ■ , 4>N-k+i) represent 
the so-called partial (or incomplete) Bell polynomials, 
which are defined as 



N-k+1 



B N>k ({<t> n })=N!j: n ^t(!t 



(13) 



where the summation takes place over all integers c„ > 0, 
such that J2 n c n = k and J^n nc « = N. 



Finally, comparing Eqs. (J6)) and (|T2j) one can readily 
obtain the following expression 

J2 / m E)g{E 1 N)e^ E dE = £ UtL Blfik ({<l> n }). 
fe=i^° fe=i 

(14) 

Given combinatorial meaning of Bell polynomials (see 
Sec. IIII A|) and also having in mind physical meaning of 
the probability f(k, E), one can find that the equality of 
the two series in Eq. (| 14|) implies equality of their terms 
(see Sec. IIII Bp . i.e. 

f(k 1 E)g{E 1 N)e-f iE dE = L^Lfl^({^}). 

(15) 

Using properties of Bell polynomials (see p. 135 in [10(), 
a k b N B Ntk ({cf> n }) = B N , k ({ab n <t> n }), (16) 
Eq. (fT5)l can be simplified to 

r°° 1 

J f{k 1 E)g{E 1 N)e-^ E dE = —B N . k {{w n }) 1 (17) 
where 

W n = ~Hn- (18) 

Eq. (dU) for the integral transform of g(E, N)e~ /3E dE, 
with the kernel function given by f(k,E), is the main 
result of this paper. The formula describes probabil- 
ity that a general system of TV particles at tempera- 
ture T — (fcs/3) -1 , regardless of its energy, consists of k 
weakly coupled clusters. Furthermore, it provides a neat 
probabilistic description of microscopic arrangements (in 
terms of clusters) of the system investigated. 

In the next section, a few points about Eq. (fT7| arc 
picked up. In particular, a novel and strictly microscopic 
interpretation of the grand thermodynamic potential is 
discussed and a new general combinatorial formula for 
the density of states is derived. The case of an ideal gas 
is considered as an example. 



III. DISCUSSION OF EQ. (fTf)) 
A. Bell polynomials 

The essential difficulty with Eq. (jTT)) arises from in- 
terpretation of Bell polynomials, B^^ k ({w n }), which are 
described by Eq. (|13[) . For this reason, we begin by elu- 
cidating their combinatorial meaning [TJ, [l7| ■ 

Suppose that N particles labelled by elements of the 
set [N] = {1,2,..., TV} are partitioned into clusters in 
such a way that each particle belongs to a unique cluster. 
Obviously, for each composition of [N] into k disjoint 
non-empty clusters (subsets, blocks) of rii > elements 
each, there are exactly 



N 

ni,n 2 , ...n k 



TV! 



n\\n2\ ■ 



■ n k \ 



(19) 
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such compositions (partitions, distributions), where n\ + 
n 2 + ■ ■ ■ + Ufe = N. Suppose further that each cluster 
of size n in such a composition can be in any one of w n 
different internal states. Then, the number of configura- 
tions of the system of N particles with k clusters is given 
by 



N-k+l 



n ffi 



(20) 



71=1 



where c„ > stands for the number of clusters of size n, 
the largest cluster size is N — k + 1 , and 



N-k+l 

E • 

n=l 



= k, 



N—k+l 

E » 



= N. 



(21) 



Eq. ([20]) describes system of N distinguishable parti- 
cles with k clusters. If one assumes that the considered 
particles are indistinguishable, and also that clusters of 
the same size are indistinguishable, the number of con- 
figurations (compositions) becomes 



N-k+l 



n c i( n i) ' 



(22) 



c ' 



which corresponds to Eq. ([20)) divided by AH 

Finally, summing Eq. (|22|) over all integers c„ > 
specified by Eqs. (|2"Tj) one gets the right-hand side of our 
main result, Eq. (|TT|) . 



B. How Eq. (Tl4)l implies Eq. ([17]) 

A pivotal role in our reasoning is played by the proba- 
bility distribution f(k, E) that the considered system of 
N particles having energy E consists of k non-interacting, 
disjoint, and statistically independent clusters. In the 
simplest case of N = 1, Eq. (fT4| is equivalent to Eq. (fT7| . 
i.e. 



poo 

/ f(l,E)g(E,l)e-P E dE = B lil ({w n })=w 1 . 
Jo 



(23) 



In words, the last expression describes probability of the 
system which consists of the only one particle. Obviously, 
in this case the particle itself is the only cluster of size 
n = 1, and its number of internal states is w\. 

In the case of N — 2, Eq. (|14|) can be written in the 
following way 

[>00 POO 

/ f(l,E)g(E,2)e-f> E dE + / /(2, E)g(E, 2)e^ E dE 
Jo Jo 



= ^a,! ({«>»}) 
= w 2 



^ 2 ({w n }) (24) 



In accordance with Eq. (|23l) . the second terms in each 
row of the above three-line expression represent thermo- 
dynamic probability of two, k = 2, disjoint and statisti- 
cally independent clusters of size n = 1, i.e. 

= ^2, 2 (M)= (J q °° f(l,E)g(E,l)e- 0E dE\ 

= / f(2,E)g(E,2)e-? E dE, (25) 
Jo 

where we have exploited a useful feature of the grand 
canonical ensemble that, for non-interacting particles, 
the grand partition function factorizes into a product of 
grand partition functions for each single particle. 

Therefore, since two particles, N = 2, may be placed 
in either two clusters of size n = 1, or a single one of 
size n = 2, the first terms in Eqs. (IM1) should stand 
for thermodynamic probability of the second realizability, 
i.e. they must describe the number of microstates of the 
system with the only cluster of size n = 2, 



f(l,E)g(E,2)e-? E dE=-B 2 , 1 ({i 



l }) = w 2 . (26) 



Similarly to what has been said about the the right-hand- 
side of Eq. (f2"3")) , it is clear the right-hand-side of the last 
equation, w 2 , is the number of internal states of a single 
cluster of size n = 2. 

Analogously, one can show that Eq. (ITT)) is true for all 
k > 3. 



C. Microscopic interpretation of the grand 
thermodynamic potential 

From the previous subsection it is clear that the right- 
hand side of Eq. (|17l) represents the number of micro- 
scopic arrangements of a system of A^ particles at temper- 
ature T which consists of k clusters of various sizes. The 
number of internal states, w n , of the cluster of size n is, 
up to a multiplicative constant, — /3, proportional to the 
nth derivative of the grand potential, </>„(/?), cf. Eq. (|18j). 

The understanding of the grand potential as the expo- 
nential generating function, 



z 

$(z,/3) = ]T— 0„(/3), 



71=1 



(27) 



for the number of internal states of n-clusters (i.e. groups 
of n interacting particles) is novel. It offers a completely 
new microscopic interpretation of this well-known ther- 
modynamic function. 

It is also interesting to note that from the microscopic 
interpretation of the grand potential, very general impli- 
cations for the stability requirements on the Landau free 
energy arise. In general, the classical stability theory (see 
Chap. 2H in Q, and Chaps. 12-14 in [2l|) places certain 
conditions on the sign of second derivatives of thermo- 
dynamic potentials of the system investigated. However, 
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given our derivations it seems necessary to expand the 
Gibbs stability requirements to higher derivatives. To be 
definite, from the non-negative definiteness of the coeffi- 
cients, w n , of Bell polynomials in Eq. (|17p . one gets the 
following conditions 



v„>i MP) < o. 



(28) 



One can speculate that the requirements may have con- 
nections with higher order phase transitions (in a sim- 
ilar way as the classical stability theory, which exploits 
second order derivatives, is used to characterize critical 
phenomena). The analogy will be further elaborated in 
our subsequent papers. 



D. The case of the classical ideal gas 

In this section, we apply our formalism to the classical 
ideal gas. In this case, the grand partition function can 
be written as, cf. Eq. (fT]), 



z 



Z G (f3,z) = l+Y, W (ZiW)) N =e 



N=l 



where 



3/2 



(29) 



(30) 



is the partition function of a single particle, A stands for 
the de Broglie wavelength, V is volume, and h is the 
Planck's constant. 

Using Eq. ® one finds that the grand potential of the 
gas is a linear function of the parameter z, 



(31) 



and the coefficients of the Bell polynomial in Eq. (|17p are 
equal 



Zi({3) and V„> 2 



0. 



(32) 



From Eqs. ([32)1 it follows that the only possible clus- 
ter decomposition of a general system of N particles de- 
scribed by the linear grand potential is the one in which 
each particle is located in a separate cluster. In fact, 
in the considered case of the ideal gas, where parti- 
cles do not interact with each other, the result is the 
only reasonable. Furthermore, since the number of in- 
ternal states of a single particle is Z\ = VX~ 3 , and 
the particles are statistically independent, one gets that 
thermodynamic probability of the system as a whole is 
Zi = Wi = BN,N{{w n })i which agrees with Eq. (|T7|) . 



E. The integral transform for the density of states 

Although in Eq. (|17l) . the kernel function, f(k,E), 
which describes the probability of the considered system 



of N particles to consist of fc weakly coupled clusters, de- 
pends on the system itself [Hj], a given explicitly /(fc, E) 
may also provide valuable information about the system. 

For example, let the kernel function be the Poisson 
distribution in the number of clusters, 



f(k,E)~P(k; (k)) 



-<fc> 



fc! 



(33) 



where k = E/q, (k) = (E)/q, and q stands for an energy 
quanta, such that [29( 



JV 



5>(fc; (fc))~^P(fc; (k)) = l. 



(34) 



k=l 



k=0 



Inserting Eq. (133)) into Eq. (fT7| one gets the following 
integral formula for the density of states 



P.P. 



g(E,N)e' 0E ;-,k 



E 
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l^B N , k ({-^n}), (35) 



where the symbolic notation P .T.J . . . ] represents the 
Poisson transform defined as 123, 01 



-f(x)dx = F(k). (36) 



P.T.[f(x);x,k] 



To elucidate how Eq. (|35|) can be used to describe en- 
ergy distribution in the system composed of N interact- 
ing particles, one should start with some remarks about 
the Poisson transform. Thus, when working with this 
transform it is important to understand, how it acts on 
arbitrary function. In some sense, it is reasonable to 
say that the transform, F(k), is a kind of weighted mov- 
ing average of the original function, fix). In compar- 
ison with the simple moving average, in which the un- 
weighted mean over some range of neighboring values is 
taken into account, in the Poisson transform the averag- 
ing is performed over the whole domain, x £ (0, oo), of 
the original function, fix), with the weights given by 
the Poisson distribution with the moving mean value 
which is equal to x. For this reason, due to proper- 
ties of the Poisson distribution, the averaging is actually 
done over the nearest neighborhood of x. The transform, 
F(k) = P.T.[f(x);x,k], looks like a smooth (or fuzzy) 
image of the original, f{x), and it is often reasonable to 
assume that [3(| F(k) ~ f(x)\ x —k- 

Therefore, simplifying, one could say that the left- 
hand-side of Eq. (|35p represents a smoothed thermody- 
namic probability, 

g(E, N)e-? E \ E=kq ~ ^B N<k {{w n }), (37) 

from which one gets the approximated general combina- 
torial formula for the density of states 



g(E,N)\ E=kq 



~Nl 



-Pjv,fe({-/3<W), 



(38) 



where q may be interpreted as the resolution parameter 
for the energy distribution. 



6 



IV. SUMMARY AND CONCLUDING 
REMARKS 

In this paper we have derived an exact combinatorial 
formula for probability that a general system of N parti- 
cles at temperature T consists of k weakly-coupled clus- 
ters of various sizes. We have also found that the grand 
potential (i.e. the Landau free energy) may be considered 
as the exponential generating function for the number of 
internal states of these clusters. Finally, we have derived 
an approximated expression for the density of states. 

We believe that the approach, when applied to sys- 
tems, such as lattice gases will provide important in- 
sights into mechanisms that lead to the occurrence of 
phase transitions. We also believe that the approach is 
the correct starting point for a future theory of higher 
order phase transitions. 
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